Sculpting microscopic magnetic flux landscapes in a E^S^CaC^Og+tj vortex lens 
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We demonstrate experimentally that the micromagnetic profile of the out-of-plane component of 
magnetic induction, B z , in the crossing lattices regime of layered superconductors can be manip- 
ulated by varying the in-plane magnetic field, Hn . Moving Josephson vortices drag/push pancake 
vortex stacks, and the magnetic profile, B z (x), can be controllably sculpted across the entire single 
crystal sample. Depending on the JT-history and temperature we can increase or decrease the flux 
density at the center and near the edges of the crystal by as much as ~ 40%, realising both "con- 
vex" and "concave" magnetic flux lenses. Our experimental results are well described by molecular 
dynamics simulations. 

PACS numbers: 74.72.Hs, 74.25. Qt 



In recent years dramatic progress has been made in 
the control of static flux structures in superconductors 
by the introduction of artificial vortex pinning sites (e.g. 
antidots and ferromagnetic dots) Q, 0, IE 0, 0. The 
next major challenge is to achieve dynamic vortex con- 
trol, so that different flux profiles can be realised in the 
same superconducting sample. Control of vortex mo- 
tion has recently been proposed and demonstrated 
0-SEJtHJIZH in so-called ratchet devices which incorpo- 
rate a spatially-asymmetric ratchet potential to achieve 
rectification of ac drives. However, a spatially asymmet- 
ric ratchet substrate is not a fundamental requirement 
for ratchet operation, and recent proposals 0, have 
described novel methods to control the motion of tiny 
particles in a binary mixture by the dragging of one com- 
ponent by the other. Here we describe the first such ex- 
perimental implementation of a flux lens in the binary 
vortex system present in highly anisotropic layered su- 
perconductors under tilted magnetic fields. 

Direct visualization 0, EE EE EE E3 has re- 
vealed that a tilted (away from the crystalline o 
axis) magnetic field penetrates the highly anisotropic 
Bi2Sr2CaCu208+,5 (BSCCO) superconductor in two in- 
terpenetrating vortex arrays, known as crossing vor- 
tex lattices l20| . One vortex sublattice consists of 
stacks of pancake vortices (PVs) aligned along the o 
axis, while the other sublattice is formed by Joseph- 
son vortices (JVs) confined between Cu0 2 layers. Su- 
perconducting currents generated by JVs deform stacks 
of PVs, resul ting in a mutual attraction between PVs 
and JVs |2( i |2 l |. This has been experimentally con- 
firmed ElEE"EE El El by the observation of PV 
chains which decorate underlying JV stacks in tilted mag- 
netic fields. JVs are usually very weakly pinned and can 
easily be driven by changing the in-plane magnetic field 



H\\ (= Hab), dragging PVs along with the m l22l . This 
joint JV-PV motion can be used to develop [l2| vortex- 
pumps, vortex-diodes and vortex- lenses which have two 
clear advantages over other existing ratchet devices, (i) 
The motion of vortices can be controlled without the need 
for nanofabricated samples with fixed spatial asymmetry, 
and (ii) the focusing efficiency can be easily varied by 
changing either the PV density (via the corresponding 
magnetic field component H z ) and/or temperature. Here 
we describe how to implement such lenses experimentally. 

Experimental results. — Our vortex lensing experi- 
ments have been performed on an as-grown BSCCO 
superconducting single crystal (T c = 91K, dimen- 
sions Immx0.75mmx50/im). The changes in B z aris- 
ing from PV lensing/antilensing were detected using a 
25/im wirewidth micro-Hall probe array patterned in a 
GaAs / AlGaAs 2D electron gas H^. The BSCCO single 
crystal was positioned directly above our sensor and se- 
cured with paraffin wax. The array has thirteen address- 
able elements, of which twelve were situated at different 
positions under the crystal and the remaining uncovered 
one acted as a reference. The sensor was driven by a 
45fj,A 32 Hz ac current and the Hall voltages detected 
with a lock-in amplifier. The out-of-plane H z and in- 
plane Hn magnetic field components are varied indepen- 
dently using a solenoid and Helmholtz coil pair (one of 
which could be precisely positioned vertically on a mi- 
crometer driven sta ge) , respectively. We use the in-plane 
'lock-in' transition [24j to align the in-plane field within 
±0.006° of the ab crystallographic planes. The exact 
alignment position is inferred from plots of at lock- 
in as a function of the height of the movable Helmholtz 
coil, in a small fixed out-of-plane magnetic field. At the 
start of lensing experiments a fixed PV density was es- 
tablished by field-cooling the BSCCO crystal from above 
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FIG. 1: (Color online) Changes in the local out-of-plane mag- 
netic induction, B z (vertically offset for clarity), versus the 
in-plane magnetic field, H\\ , measured near the center (a) and 
near the edge (b) of the sample at 77 K. Arrows indicate 
the direction of the Hu field sweep. The left hand inset in (b) 
shows sketches of PV density profiles across the crystal for in- 
creasing Hu (bottom to top). The right hand inset in (b) plots 

the vortex lensing amplitude (AB Z = BT ax (H l} )-BT in (H\\)) 
at H z —5.3 Oe for the central sensor as a function of temper- 
ature. 



T c in a known value of the out-of-plane field, H z . The 
in-plane magnetic field Hu was then cycled several times 
until a steady-state loop was obtained. During each cy- 
cle, H\\ was slowly (1.7 Oe s _1 ) ramped up to a max- 
imum of 150 Oe and down to a minimum of -150 Oe 
and back to zero, while the Hall voltage was monitored 
at a chosen element to measure the magnetic induction, 
B z . In practice the observed lensing was a function of 
the measurement position under the crystal. For the 
sake of brevity we only present data here for two ele- 
ments, one at the sample center and one 225/xm from 
one of the edges parallel to H\\, which fully illustrate 
the extremes of behaviour seen. Fig. ^ shows lens- 
ing data measured (a) at the central location and (b) 
near the edge of the sample at 77K for different values 
of H z . In all cases the data have been symmetrized 
(B z T (H ]{ ) = ±(B Z t {H\\) + B z I (-iJ, t ))) to ac- 
count for a small misalignment between the plane of the 
BSCCO crystal and the Hall probe array. In both cases 
for H z < 2 Oe, the PV system shows a weak reversible 
response which inverts when H z is inverted, attributable 
to the dragging of PV stacks which are all trapped on JV 
chains. For H z > 2 Oe, free PVs exist between chains 
(mixed chains/lattice state) and we start to see stronger, 
irreversible behaviour related to the compression of free 
PVs and their cutting through JV stacks at high in-plane 
fields. At higher H z the magnetization loop B z (Hu) for 
the central element has a "butterfly" shape exhibiting: 
(i) a fast increase of PV density when Hu increases from 
zero, followed by weaker (saturation-like) dependence of 
Bf ntral (if||); (ii) a rapid reduction of Bf ntral (iJ|| ) when 
H\\ decreases from its maximum value, followed by a re- 
markable antilensing (an overshoot in the reduction of 
PV density) effect Bf ntlaX (H\ { > 0) < Bf ntral (iJ|| = 0). 
Note that on the Hu -increasing branch of the lensing 
loop, Bf ntr&1 (H\\ ) is always higher than on the decreasing 
branch on the right-side (Hu > 0) of the hysteresis cycle, 

i.e., Bf nUal (H h dH\\/dt > 0) > Bf ntral (i?||, dH\\/dt < 
0). We denote such loops as "clockwise". 

The data from the Hall element near the sample edge 
(Fig. \T]p) provide insights into the spatial distribution of 
the PV density in our lensing experiments. In stark con- 
trast to Fig. we now see strong antilensing behaviour 
for H z < 6.5 Oe. This is easily understood in terms of the 
PV profiles generated during our experiments. Since PVs 
are pushed from two opposite edges of the sample towards 
the center, there must be regions near these edges which 
experience a decrease in PV density, while accumulation 
is occurring in the crystal center (see sketched profiles in 
the left hand inset of Fig. [TJj). It is interesting to note 
that the counter-clockwise loops at low H z transform to 
clockwise ones when H z increases. This is best illustrated 
in the curve at H z = 6.5 Oe, which shows an extra cross- 
ing point for ify > between traces on the sweep-up and 
sweep-down. At such high values of H z PVs penetrat- 
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ing through the sample surface partially compensate for 
the deficit of PVs near the edges, and enhanced PV-PV 
repulsion results in a broadening of the focus region and 
a shift of the lensing/antilensing interface towards the 
sample edge. Hence a transition from counter-clockwise 
to clockwise B z loops occurs as the out-of-plane field is 
increased. We conclude that the PV profiles, and their 
changes, are a subtle function of the applied out-of-plane 
field, allowing us to precisely control the "magnetic land- 
scape" inside the sample. Vortex control is achieved here 
using just the so-called dc-driven (quasi-adiabatic) mode; 
ac-drives with trains of time-asymmetric Hn pulses will 
be described elsewhere |2Sj. 

The right hand insert of Fig. ^ shows a plot of the 
lensing amplitude at H z =5.3 Oe as a function of temper- 
ature. Surprisingly, we find that the lensing amplitude 
falls rapidly with temperature, being very weak at 84 K 
and undetectable at 86 K. Since this is well below the 
critical temperature of the BSCCO crystal (T c ~ 91K), 
and the crossing lattices interaction strength is only very 
weakly temperature dependent |2flj . we conclude that 
efficient lensing requires finite bulk pinning to prevent 
PVs escaping from focus regions laterally parallel to JV 
stacks. The effectiveness of bulk pinning will drop rapidly 
at elevated temperatures leading to a rapid reduction in 
efficiency. 

Simulations. — The minimal model to simulate the ob- 
served lensing effects describes the overdamped dynam- 
ics of JV and PV rows using a set of coupled equa- 
tions of motion: ^rjjxf /a J = ff J + ff H + ff p , and 
rj P x P /a p = f PP + f PH + f PJ , where xj and x P are 
the positions of JV and PV rows with distances between 
JVs(PVs) in a row a J /"f(a p ) and anisotropy parameter 7, 
while rjj and rjp are the JV and PV viscosities. The vis- 
cous forces slowing down the vortex motion are balanced 
by: (1) the repulsive force f JJ between JV rows (includ- 
ing images of rows with respect to the sample surface); 
(2) the interaction f JH of JV rows with Meissner cur- 
rents generated by the externally-applied time-dependent 
magnetic field H\\(t); (3) the repulsion f pp between PV 
rows (including images); (4) the interaction f PH of PV 
rows with H z ; and (5) the attractive forces / between 
rows of JVs and PVs. 

The interaction between vortex rows decays expo- 
nentially [il: ff J r/Dr]j = (a p '/3 / a J ) J2j sgn[a;/ - 
x/]exp[-|a;/ - xj\/X c ], f^r/Dr/p = £\ sgn[x P - 
x p ] exp[— \x P — x p \/X a b], where the interaction lengths 
are the in-plane A a b and c-axis A c penetration lengths, 
while (3 — rjpK'Yqj) is the relative viscosity. Hereafter, 
we normalize all distances by the half-width of the 
sample D and all time scales by r = lQnr\ 2 b a p r]pD /$q. 
The distances between PV and JV rows are related 
to the H z and Hu magnetic field components by 
a p « (Qq/H,) 1 / 2 and a J fa (7$ /-H|| ) 1/2 - The interac- 
tion with the Meissner current decays exponentially on 
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FIG. 2: (Color online) Simulated loops of the local induction 
B z {Hu) at the center (a) and near the sample edges (b) for 
different values of H z (H z increasing from bottom to top in 
each panel). The red region in each inset indicates where B z 
was monitored in the sample. The main features found in 
experiments (Figs. ^,0d) are reproduced in the simulations. 



the scales A c and X a b from the surface (x = ±1) of the 
sample, for JVs and PVs, respectively. The correspond- 
ing forces can be modelled (2{J as fi JH T/Dr]j = 
-(2A c a p /3/(a J ) 2 )sinh(a:,f/A c )/cosh(L)/A c ) and 
f PH T/D VP = -{2\ ab /a p )sinh{x P /X ab )/ cosh(D/\ ab ). 
The JV-PV attraction can be approximated 
as ff p r/Dnj = X 2 a M(a J ) 2 £ fe sgn(^ 
xf ) exp(-47r|a;f - xj\/a J ) and f PJ r/Dr]p = 

(A> p /? 1 /(a^ 3 )E I sgn(x/- a; nexp(-47r|^-^|/ a J ), 
where A w 167r 2 /ln{l + [{X 2 Ja J ) 2 + l]/[(A ab /a p ) 2 + 1]} 
is related to the tilt elasticity of PV stacks. The 
qualitative agreement between the very complex and 
nontrivial experimental data and simulations indicates 
that the model used captures the essential physics. 

Our molecular dynamics simulations follow the experi- 
ments. First H\\ was slowly increased from zero to i?| nax , 
and then decreased back to zero over the same period of 
time. Such cycles were repeated several times to achieve 
steady-state loops. The average PV density was moni- 
tored (as a function of Hh) both at the center and near 
the edges of the sample. 

Comparison theory- experiment. — The simulated 
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FIG. 3: Net change of magnetic induction, AB Z , as a function 
of H z measured at the sample center for four different tem- 
peratures. Inset shows the simulated lensing efficiency. Both 
experiment and simulation exhibit a maximum as a function 
of H z . 



clockwise B z (Hn) "butterfly" loops for the PV density 
at the center of the sample show the same features |27j 
that were observed in experiments (c.f. Figs. ^ and|30 
and can be easily interpreted: (i) First, B z entTal increases 
with H\\ as JVs move towards the sample center and 
drag PVs with them. This is consistent with theoretical 
predictions (ii) At a certain in-plane field H», the 
PV density at the center of the sample saturates and 
even starts to decrease. The PV density at the center 
is now large enough that PV-PV repulsion becomes 
very large. PVs start to cut through the JV rows and, 
near the maximum compression point, some of these 
PVs escape the narrow potential produced by the JVs 
and are lost; (iii) On the decreasing branch of the loop, 
both experiments and simulations exhibit a remarkable 
anti- lensing effect. This arises because a smaller total 
number of PVs now spreads out over the whole sample, 
resulting in a deficit of PVs at the center. Also the ratio 
of the lensing to antilcnsing effect, max[£^ cntlal (-ff||) — 
Sf ntral (0)]/|min[Bf ntral (-ff||) - £^ cntral (0)]| decreases 
to about one when the out-of-plane field increases, in 
agreement with experiments (this produces rounder 
loops for higher H z , see Figs. ^,|^,). 

The simulations also capture all the main features 
found in the experiments near the sample edge: (i) At 
very low out-of-plane fields, H z , the area of this re- 
markable counterclockwise loop increases with H z ; (ii) at 
higher out-of-plane fields the counterclockwise loops nar- 
row, transform into clockwise loops and broaden again. 

Fig. |2| illustrates the experimental (main figure) and 
simulated (inset) dependence of lensing amplitude on 
out-of-plane magnetic field. In both cases the amplitude 
shows a pronounced peak as a function of H z . At 77K 



the experimental lensing efficiency exhibits a maximum 
of nearly 40% at H z ~ 5 Oe. This field value is in rea- 
sonable agreement with the predicted maximum in JV 
pinning strength (by PVs) at B z ~ 0.26^, ~ 6G |2s|. 

Conclusions. — We have experimentally realized con- 
cave and convex magnetic vortex lenses which fo- 
cus and defocus the out-of-plane magnetic flux in a 
P^S^CaC^Os+a sample. Remarkably, this was done 
by employing the "dragging" of PVs by JVs with no fixed 
spatial asymmetry in the sample. We show that the PV 
density near the center/edges of the sample can easily 
be controlled by changing either the in-plane or the out- 
of-plane fields, as well as by varying the temperature. 
The experimental results are well described by a simple 
model considering the dragging of one vortex species by 
the other. This novel method of quantum-motion-control 
(a ratchet with no ratchet potential) opens up new av- 
enues for the manipulation of flux quanta and nanoscale 
particles. 

We acknowledge support from EPSRC in the UK un- 
der grant No. GR/R46489/01, the ESF VORTEX net- 
work, the US NSA and ARDA under AFOSR contract 
No. F49620-02-1-0334, NSF grant No. EIA-0130383, and 
a Grant-in- Aid from MEXT, Japan. 



[1] M. Baert et al, Phys. Rev. Lett. 74, 3269 (1995). 
[2] J.I. Martin et al, Phys. Rev. Lett. 79, 1929 (1997). 
[3] D.J. Morgan et al, Phys. Rev. Lett. 80, 3614 (1998). 
[4] M.J. Van Bael et al, Phys. Rev. B 59, 14674 (1999). 
[5] I.F. Lyuksyutov et al, Phys. Rev. Lett. 81, 2344 (1998). 
[6] J.F. Wambaugh et al, Phys. Rev. Lett. 83, 5106 (1999); 

C. J. Olson et al, ibid 87, 177002 (2001); B.Y. Zhu et 

al, ibid 92, 180602 (2004). 
[7] W.K. Kwok, et al, Physica C 382, 137 (2002). 
[8] J.E. Villegas et al, Science 302, 1188 (2003). 
[9] J. Van de Vondel et al, Phys. Rev. Lett. 94, 057003 

(2005). 

[10] Y. Togawa et al, Phys. Rev. Lett. 95, 087002 (2005). 
[11] R. Wordenweber et al, Phys. Rev. B 69, 184504 (2004). 
[12] S. Savel'ev and F. Nori, Nature Materials 1, 179 (2002). 
[13] S. Savel'ev et al, Phys. Rev. Lett. 92, 160602 (2004). 
[14] T. Matsuda, et al, Science 294, 2136 (2001). 
[15] A. Grigorenko et al, Nature 414, 728 (2001). 
[16] V.K. Vlasko-Vlasov et al, Phys. Rev. B 66, 014523 
(2002). 

[17] M. Tokunaga et al, Phys. Rev. B 66, 060507(R) (2002). 
[18] C.A. Bolle et al., Phys. Rev. Lett. 66, 112-115 (1991). 
[19] L.N. Bulaevskii et al, Phys. Rev. B 46, 366 (1992). 
[20] A.E. Koshelev, Phys. Rev. Lett. 83, 187 (1999). 
[21] S.E. Savel'ev et al, Phys. Rev. B 64, 094521 (2001). 
[22] G.K. Perkins et al, Supercond. Sci. Technol. 18, 1290 
(2005) 

[23] M.S. James et al, Phys. Rev. B 56, R5771 (1997). 
[24] A.E. Koshelev, Phys. Rev. B 48, 1180 (1993). 
[25] D. Cole et al, Nature Materials (in press). 
[26] S.E. Savel'ev, V.S. Gorbachev et al, JETP 83, 570 
(1996). 



[27] Animations illustrating vortex lensing are available at [28] A.E. Koshelev, Phys. Rev. B 68, 094520 (2003). 
http:/ /dml.riken.go.jp/vortex-dc 



